CONSTRUCTION OF EXACT TRAVELLING WAVES FOR THE 
BENJAMIN-BONA-MAHONY EQUATION ON NETWORKS 
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Abstract. We are interested in the existence of travelling waves for the Benjamin- 
Bona-Mahony equation on a network. First we construct an explicit wave, de- 
fined in M. Then, we use this wave to derive some conditions on the coefficients 
appearing in the equations and on the geometry of the network to ensure the 
existence of travelling waves on the network. 



1. Introduction 

Partial differential equations on networks (also known as metric graphs or quan- 
tum graphs) have been studied for many years now. While they originate from the 
pioneering studies of quantum chemists Klaus Ruedenberg and Charles W. Scherr 
in [5D] , the current state of the art of this theory owns much to the thorough math- 
ematical analysis performed on these problems in the 1980s in particular by Giinter 
Lumer and some of his students and collaborators. These and subsequent mathe- 
matical investigations have then paved the road to the re-discovery of this kind of 
models in the context of theoretical quantum physics in the last decade. Perhaps 
because of the quantum physical bias in the most recent investigations, little atten- 
tion seems to have been devoted to the study of nonlinear phenomena in networks 
- and anyway, mostly from a linear point of view. However, in the theory of non- 
linear PDEs one is often interested in questions that are somewhat complementary 
to those considered in the linear case. For instance, one is often particularly keen 
on finding out whether a given system modeled by a nonlinear PDE can support 
any travelling wave, for certain conveniently chosen initial data (which then turn 
out to be suitable inputs), rather than asking for mere existence and uniqueness 
of a (possibly non-physical) solution for any initial data. In the present paper we 
are going to suggest a possible general approach to the study of travelling waves in 
networks. It turns out that, unlike on the case of PDEs defined in the whole space 
(but somewhat similarly to the case of PDEs on Riemannian manifolds 0) exis- 
tence of travelling waves imposes severe restrictions on the geometry of the graph. 
The ideas presented here can be applied to a large class of nonlinear equations 
featuring a second order linear differential operator as the leading order term, in- 
cluding the nonlinear Schrodinger equation [T], the one-dimensional Navier-Stokes 
equation [12], or the FitzHugh-Nagumo or Rail equations [THl [TO], as soon as a 
special solution of the PDE without boundary conditions is known. For example, 
we may as well borrow the travelling wave analysis from |13j in order to discuss the 
porous medium reaction-diffusion equation on networks. 
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In this note we prefer to fix the ideas and focus primarily on the Benjamin-Bona- 
Mahony equation 

(BBM) Ut — aur^j-t + buuj. + duj. = in M for t > 0, 

where a € (0, oo), 6 g M* and d G M. This equation models long waves in non-linear 
dispersive systems and is known as a good substitute for the Korteweg-de Vries 
equation 

Ut — Uxxx + =0 in M for t > Q 

in the case of shallow waters in a channel, see [22 . In [6 , Benjamin, Bona and 
Mahony studied the initial value problem corresponding to (BBM) and they estab- 
lished global existence and uniqueness results. 

In [8^ , Bona and Cascaval considered a finite metric tree consisting of edges and 
established the well-posedness of a system of BBM equations 

(1.1) Ut — aiUxxt + biUUx + diUx =0 on each for t > 0, 

with standard continuity conditions and Kirchhoff conditions at vertices v = H &j 
and with Dirichlet boundary conditions at endpoints. This equation is used to 
model the blood flow in the human cardiovascular system, see ^\ and references 
therein. However, just thinking of the cardiovascular system seems to motivate the 
discussion of topologies that include circuits. 

We begin our note discussing a general way of deriving conditions on the coefficients 
of (BBM) as well as on the metric and orientation properties of the network (but 
not on the topological ones!) as a consequence of the standard boundary conditions 
that are customarily imposed on evolution equations on networks. In Section [3] we 
show that this already excludes certain network configurations. This analysis is not 
really restricted to any specific semilinear PDE of second order (in space). 
As soon as one focuses on the BBM equation, however, it seems natural to expect 
that the pressure profile follows a wave-type behavior. Thus, our interest in this 
note is to derive an explicit formula for the travelling waves solution of (BBM) 
extending the bifurcation method (Section |4]) used by Song and Yang in the case 
of the Zakharov-Kutnetsov-BBM equations 21j. Then, in the case where the BBM 
equation is posed on a network, we use this formula to determine some conditions 
on the coefficients appearing in the equations and on the geometry of the network 
to ensure the existence of travelling waves on the network. More precisely, we 
construct a solitary wave u such that u(t,x) — ip{x — ct) on each edge with some 
if vanishing at ±oo. For the reader's convenience, we begin with the case where 
the network is a star (Section [s]), then we deal with the case of a tree (Section |6]), 
and finally we treat the general case of a network that may possibly contain circuits 
(Section [f]). 

It seems that not many other investigations on travelling waves or solitons on net- 
works have been performed. In fact, we are only aware of [H Sections 16-17] and [S] 
for the case of linear diffusion and a certain class of reaction-diffusion equations; and 
some recent progresses in the theory of linear and nonlinear Schrodingcr equations 
on star graphs, see [Tj and references therein. 

2. Preliminaries on networks 

For any graph F — (V, E, e), the vertex set is denoted hy V = V{r), the edge set 
hy E — E(r) and the incidence relation by GC V x E. All graphs considered in this 
paper are assumed to be non-empty, simple, connected and at most countable. The 
simplicity property means that F contains no loops, and at most one edge can join 
two vertices in F. We give a numbering of the vertices v.^, i G V C N with card V = 
card V, and a numbering of the edges e^, j G E C N with card E = card E. We 
denote by 
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7V(v) := {j e N / V e e,-}. 

the of boundary index set of v, i.e., the set of indices of those edges incident in v; 
while d(v) denotes the degree of v, i.e., 

d{y) := card N{\/) 

Also, we assume our graphs to be locally finite, i.e., 

V V e V(r), d{y) < 00. 



Recall that a simple, connected graph is called a tree if it does not contains any 
circuits, i.e., if any two vertices arc connected by exactly one undirected path; it is 
called a star if it is a tree and all vertices but one have degree 1. 

We then consider networks by associating with a given graph a topological struc- 
ture in M™, i.e., we regard V^(r) as a subset of M"* (in fact, it is well-known that each 
countable graph can be embedded in M^) and the edge set consists of a collection 
of Jordan curves 

EiT) = {7rj:[0,ej]^M"' / je-E} 

with the following properties: each support Sj := ttj {[0,£j]) has its endpoints in 
the set V"(r), any two vertices in ^(r) can be connected by a path with arcs in 
E{T), and any two edges e^- / eh satisfy Sj n eh C V^(r) and card(ej- fl eh) < 1. The 
arc length parameter of an edge is denoted by xj. Unless otherwise stated, we 
identify the graph T = (V, E, e) with its associated network 

G=\j7rj{[0,£j]), 

especially each edge nj with its support e^-. G is called a C-network, if all tTj G 
C'^{[0,£j],M."^). We shall distinguish the boundary vertices 

Vb = {vi€V / d{vi) = 1} 

from the ramification vertices 

Vr = {vieV / d{\Ji) > 2}. 

The orientation the network is given by means of the incidence matrix I := (oij), 
where 

{1 if 7rj(£j) =Vi, 
-1 if^,(0)=Vi, 
otherwise. 

The cases of t^j = 1 and ty- = — 1 correspond to the cases of an incoming and an 
outgoing edge, respectively. 

The reason why we speak of networks, at the risk of confusing a reader more fa- 
miliar with graph theory, is that our setting is in fact a generalization of the graph 
theoretical theory of networks, where the edge lengths can be seen as their capaci- 
ties. Accordingly, wc also borrow further graph thcioretical notions and speak of a 
sink (resp., a source) to describe a vertex with only incoming (resp., only outgoing) 
edges. 

For a function m : G x M ^ K+ we set Uj{-,t) := u{Trj{-),t) : [0, ij] M for alH > 
and use the abbreviations 



d 

Uj{t,Vi) := Uj{TrJ^{\/i),t), djUj{t,\/i) := -K—Uj{t,Xj) 



dxj TT-^K) 



etc. 
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3. The BBM equation on networks 

Our analysis is based on the classical observation that, by definition, existence 
of a travelling wave solution for a separable evolution equation 

(3.1) F{t,x,v)^0, 

(where F may possibly depend on partial derivative of any order oiv), i.e., existence 
of some function -0 and some constant c > for which the solution of (3.1 1 satisfies 

(3.2) v{t,x) = ip{x - ct), 
is equivalent to solvability of the system 

F{t,x,v) = 0, 



(3.3) 



because by (3.2 1 the chain rule applied to ^jj yield 

Uf (i, x) — — C^'(x — Ct) — —CVx(t, x), 



and conversely ( 3.2 ) yields the only possible solutions of the second equation of ( 3.3 1 . 
In particular, in the case of the (BBM), this Ansatz leads to considering the system 

, , ( Ut - aUxxt + buUx + dux = in M for i > 0, 

^ ' ^ \ Ut = -cux in M for t > 0. 

Solving this system amounts to finding a function w of one variable along with some 
wave velocity c £ K such that 

(3.5) acip'" + bipLp' + {d ~ c)ip' = 

where (p' denotes the derivative of (p, with a > 0, 6 7^ and d e M. If such a (p 
exists, then by definition u will be obtained by 

u(t, x) :— if{x — ct), 

on each edge. Still, existence of a solution of ( |3.5| is a necessary but not sufficient 
condition for the whole network to support a travelling wave. In order to construct a 
travelling wave solution one needs to transform boundary conditions for m, which we 
will introduce soon, into conditions for p. We will do so applying an idea developed 
and thoroughly discussed in [4, Sections 16-17]. 

To fix the ideas, as the basic geometric transition condition at ramification ver- 
tices we impose the continuity condition 

(3.6) Vv, e Vr, Vi e M+ : e^- n e,, = {vj =^ Uj{t,y{) = u,,(i, v,), 

that clearly is contained in the condition u G C(G' x M). Moreover, at all vertices 
Vi e Vr wc impose the classical Kirchhoff condition 

(3.7) Vvi e Vr, Vt e M+ : ^ iijajdjUj{t,v,) = 0, 

where aj is the coefficient of the BBM ( |1.1[ ) on the edge . This Kirchhoff conditions 
corresponds to imposing conservation of the flow - and hence of the mass - at each 
ramification vertex. Including the coefficients in this condition is necessary in 
order to make this conservation independent of the length of the edges (which in 
turn depend on their parametrization) . Note that Condition (3.7) does not depend 
on the orientation. 

Summing up, in the present paper we consider a system of BBM equations on a 

C^-network G 

(BBMG) 

dtUi - aidfdtUi + biuAui + didiu^ = 0, x e e^, t > 0, 

Uj{t,Vp) = Uk{t,\/p) j, k e N{\/p), t > 0, 

LpjajdjUj{t,Vp) =0 t >0, 
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where a.j > 0, 6^ e M*, G M, for alH e E and all Vp e 

Initial data are not prescribed, since such data would already fix the initial profile 
of the front wave. Also, we do not impose boundary conditions on the boundary 
vertices v € Vf,, since such data describe the tail of the front wave. 

Remark 1. Another approach would be to replace each edge e containing a vertex 
in Vb by a half-line whose endpoint is enVr- In this way we would consider a kind 
of nonlinear scattering problem. 

Definition 3.1. A strong solution of the system (BBMG) is a function 

u e I w e C(G) I \/ ieY, Ui£ C^'^(e, x [0, oo)) and dtu, e C^'^ie^ x [0, cx)))| 

that satisfies (BBMG) pointwise. 

More specifically, in this paper we are looking for travelling wave solutions. 

Definition 3.2. A strong solution u o/(BBMG) is called a travelling wave if there 
exists a velocity vector (ci)igE C and a vector-valued function {ipi)i£E C C'^(M) 
such that 

Ui{xi, t) ~ ipi{xi — Cit) for all Xi G and t>Q. 

Definition 3.3. A travelling wave u defined by 

Ui{xi,t) :— ipi{xi — Cit) for all Xi G and t >0, 

is said to be stationary if on each edge e,j either ip'^ = or Ci — 0. We call ip solitary 
if it admits at most one local extremum and if lim ip{z) exists in M. 

Observe that the development of a travelling wave in a network on each of whose 
vertices a boundary condition is imposed requires the existence of paths of infinite 
length, possibly allowing repetition of edges but not doubling back. Hence, it is 
apparent that only infinite graphs and/or graphs with circuits can support a travel- 
ling wave - unless we drop any condition on the function at the boundary vertices. 
It turns out that in fact additional compatibility conditions are necessary. 

Such conditions can be derived by the standard boundary conditions (continuity 
conditions (3.6) and Kirchhoff conditions (3.7|) which we impose at the ramification 
vertices. This idea has been exploited already in [4]. Indeed, if we already assume 
the solution to be a travelling wave, then the continuity condition (3.6) is equivalent 
to 

(3.8) Vv,; e K Vt € K"*" : e^- n = {vj =^ (pj{eij - cjt) = (pk{eik ~ Ckt), 
while derivating both members of (3.6) with respect to time we see that 

(3.9) Vvi e KVi e : e^- n = {vj ^ Cj'^'ji^i] - Cjt) = Ck^p'ki^ik - Ckt) , 
with 



Ml±^. if..(^.) = v., 

2 |0 if^,(0)=v„ ''^ 



Remark 2. Using the variable z = x — Cjt we see that (3.8) can be equivalently 
re-written as 

(3.10) Vvj e Vr : ej Ciek = {vj => ipj{z) = ipk\£tk - + —zj. 

\ Cj Cj I 

Hence, the continuity at ramification vertices implies that for any fixed pair of mu- 
tually adjacent edges Sj , e^; , each ipj is of the form 

(3.11) (p,(z) = ^fe(c(fc,j) + %). 
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for a constant C'{k,j). But then, owing to connectedness of the graph, (3.11) can be 



extended to any pair of edges &j,e.k, where C{k,j) is a constant depending on the 
edge lengths and speeds along a suitable path containing e.j and . This shows that 
a travelling wave is completely determined by its profile on one single edge and 
by the speeds Ci, C2, . . . Using (3.10), we finally observe that 

Ck 



n Gfe {vj djU.j{t,Mi) = — (^'j,(ejfe - ckt), 



t > 0. 



Combining this relation with (3.9) we see that a travelling wave u satisfies the 



Kirchhoff condition (3.7) at a ramification vertex Vj if and only if 

a. 



(3.12) 



= 0. 



We remark that, unlike (3.9| and (3.8), the compatibility condition (3.12) does not 



impose any restriction on the geometry of the network, hut only on the coefficients 
appearing in (BBM) on adjacent edges. We also remark that the system of ODEs 
consisting of the equation (3.5) on each edge is underdetermined, since the leading 
term ip'" is of third order but we are only imposing conditions (3.8) and (3.9) on if 
and Lp' . 

Lemma 3.4. // a travelling wave solution of (BBMG) is stationary, then it is 
constant. 



Proof. By Definition 3.3 the travelling wave u is stationary if and only if on each 
edge Sj either Lpj is a constant function or Cj vanishes. In either case, Lpj will be 
constant in time (see (3.5 1 ) on each edge e.j. In view of the continuity conditions 
(3.6) and because G is connected by assumption, the solution u will be constant in 
time, as well. □ 

Thus, we focus on the case where the wave u really travels, i.e. u is non-constant 
and all the Ci are strictly positive. 

Lemma 3.5. If a non-constant travelling wave solution o/(BBMG) exists, then no 
ramification vertex can be either a sink or a source. 

Proof. Since we restrict ourselves to the case Oj > and Cj > 0, if a non-constant 
travelling wave exists in G, then it follows from (3.12) that neither 

a," 



> 



(as it is the case, in particular, for a sink - i.e., if i.y > for all j E E) nor 



J6E 



< 



(as, in particular, for a source - i.e., if t^j < for all j G E ) can occur at any vertex 
Vi e Vr- □ 

We have actually proved slightly more: If in fact aj — cj for all j £ E, then 
the above lemma states that each vertex has to be balanced, in the sense that each 
vertex has to have the same number of outgoing and incoming edges. This condition 
also appears in the theory of first order differential operators on network [17 , and is 
known to be necessary for the existence of travelling wave solutions for a nonlinear 
Schrodinger equation on a star graph p . Of course, in classical graph theory it is 
well-known that a directed graph is Eulerian if and only if it is balanced. 

Remark 3. Just as continuity and Kirchhoff-type conditions are natural for a man- 
ifold of partial differential equations on networks, and not only for the (BBMG) 
system, the compatibility conditions (3.8) and (3.12) are the natural one for all 
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problems on networks involving differential operators in divergence form. Results 
similar to those of this note could then be deduced for different classes of evolution 
equations. 

In fact, even if we have chosen to concentrate on the BBM equation, the above 
analysis shows how our strategy can be applied to more general problems, as those 
mentioned in the introduction. Though, other favorite differential models for waveg- 
uides - like the KdV equation, the Camassa-Holm or the Whitham equation - do 
not seem to be treatable by our methods, as it is not quite clear which boundary 
conditions should be imposed on a third order differential equations or on integral 
equations [9l[2l[T4]. 

4. Profile of the front 
We want to determine an explicit function : M — !■ M such that u defined by 
u{t,x):^ip{x-ct), (i,x)eMxM+, 
for some c > 0, is a solution of (BBM). Integrating ( |3.5[ ) over x, we are lead to 

—c^p + aap" + ^f"^ + dtp = A, 

for some A e R. With the notation ip (p' , we obtain the first order differential 
system 



(4.1) 



(p' = ip, 



ac 



As in we study the phase portraits of this system. First, we consider the 
functional 

(4.2) Hi^, ^) ^ ^ _ 1 + _ b \ 

2 ac\ 2 / 



A direct calculation shows that H is constant along any trajectory of (4.1), i.e. for 
all solutions y i— >■ {tp{y), 4'{y)) of system ( |4.1[ ) we have 

j^H{^{y),i^{y))=0. 



Then, we need to investigate the stationary points of (4.1). 
• If 

(c - df + 2Ab < 0, 



then (4.1) has at most one stationary point and all its trajectories are un- 
bounded, see Figure [T] 

Hence, we focus on the case (c — rf)^ + 2Ab > 0. 

• If 

{c-df + 2Ab>0 
then (4.1) admits two stationary points 



fc~d- J{c- df + 2Ab 
P. = [ ^ ,0 

and 

'c- d+ Jic- dy + 2Ab 

P2- ' 



Clearly, the eigenvalues A of the linearized system of (4.1 1 around pi satisfy 



2 _ V(c-d)2 + 2A6 
A — 
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6 > and (c - d)^ + 2Ab^0 



6 < and (c - d)^ + 2Ab = 



Figure 1. Phase portraits of (4.1 1 with (c - dy + 2Ab < 0. 



and the eigenvalues /i of the Hnearized system of (4.1 ) around pi satisfy 

y/{c- d)^ + 2Ab 



ac 



Then according to the theory of dynamical systems (e.g. [3], [H] and [15]). 
we obtain that pi is a saddle point for (4.1), whereas P2 is a center. 



.3) 
•4) 



Again in the case 

(c - df + 2Ab > 0, 
we are specially interested in the homoclinic orbit F. The heteroclinic ones 
also represent travelling waves, but we can not derive their explicit formula 
(see Remark [5|, and some of them are singular at ±oo (e.g. the branches 
El and E2.). Indeed, a homoclinic orbit corresponds to a solution (<^, "0) of 
the (4.1 1 defined on M and satisfying 

lim f{y) = - and lim ^{y) — 0. 

y— f±oo b y— ^±00 



In the special case A — 0, (4.3) reduces to 

lim fiv) = lim V'(y) = 0, 




(4.5) 



and recalling that H is constant along any trajectory of (4.1 1, we obtain 



H(if, i/)) = along F, and in the ip — ip plane, using (4.2 ), T can be described 
as 

c — d ry b , 



r 



ac ' 3ac 

Up to a translation, we can suppose that ((/9(0), ■0(0)) = (iy9o,0). Then, 
from standard regularity results (see [Hj) the abscissa of trajectory T 
is a solution of (3.5) belonging to C°°(IR). When 6 > 0, (yS is positive in 
M, increasing in (— oo,0) because i^{t) — (p'{t) lies in upper half-plane for 
t < 0, and decreasing in (0, oo) because lies in lower half-plane when t > 0. 
When 6 < 0, <p is negative in K, decreasing in (— oo, 0) and increasing in 
(0,oo). In both cases, it is easy to see that </? is a solitary wave of (3.5), see 
Figure [3] 





Figure 3. Solitary wave of (3.5 1 



From (4.5 1 and writing ip — (p' — 'j^, we obtain 



ac > Sac > 



dp = dy, 



and integrating between and y, we have 

sgn(^)A/-^^ / / ^ . ^=ds = / dy = -y, 



c — d 



viv) Ws 



_b „3 



3(c-d)" 



for y < 0, 



V 
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and 



-sgn(6) 



v{y) 



-ds 



dy = y, 



for y > 0. 



In order to compute explicit formulas, we need to impose 
(4.6) Od. 

Using both anti-derivatives x H> arccosh^^ — ij and x i~> arccosh^^ + 1^ with 



a 



± 



e(c-d) 



and up to a translation, we obtain: 



Theorem 4.1. Let A = 0. When a > 0, b ^ 0, d e 
admits the solitary waves 

6{c-d) 1 



(4.7) 

and 
(4.8) 



viy) 



f{y) 



1 + cosh 



6(c - d) 



1 — cosh 



and c > maxjd, 0}, (3.5) 
/or y e M 

foryc, M*. 



Remark 4. The smooth solitary wave Lp describes the homoclinic orbit T, and the 
singular solitary wave ip describes the branches Si and S2- 



If we do not impose ^ = in (4.4), we obtain a solitary ipA wave solution of (3.5) 
satisfying (4.3). Let (^q = — A. From, (3.5), we obtain that 

(4.9) 



c(Pq + acip'^' + b(po(pQ + {d + A)(p'q = 0, 



and ifiQ verifies (4.4). Hence, from Theorem 4.1 we have: 



Corollary 4.2. Let A e 



d- A)2 4 

VA{y) 



If a > 0, b ^ 0, d eR and c > max{(i 
4A6 > 0, (3.5) admits the the solitary wave solutions 

&{c-d + A) 1 



A 



1 + cosh 



c-d+A 



A, 0} satisfy 
for y eR 



iPAiy) = A + 



6(0- d + A) 



1 — cosh 



for y e 



Remark 5. On can see in Figure [1| that there also exist some periodic orbit of 
(4.1). It corresponds to a periodic wave for (3.5), but we are unable to derive 
explicit formulas, even if it is possible to express it in terms of the Weierstrafl's 
function, see [El [23] . 

5. BBM EQUATION ON A STAR 

Let us consider a semi-infinite star, i.e., the finite graph with one vertex Vi and 
N edges of semi-infinite length. In view of ( |3.12 ), we suppose that the incidence 
vector (tii)i<i<Ar, defined in (2.1), is not ±(l)i<i<Ar. Thus, up to relabeling, there 
exists 1 < L < iV such that 



1 for l<j<L, 

-1 for L + l < j < N. 

In other words, we have L incoming edges ei, . . . ,6^ and N — L outgoing edges 
e^-i-i, . . . , ejv- We identify the incoming edge e^- with the half-line {—ooj, 0], and the 
outgoing edges efe with the half- line [O,oofe). 
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Figure 4. A star with infinite edges. 



We want to construct a solution u in the form Ui{t, x) = ipi(x — Cit) on each ed 
(5.1) 



ei, where Lpi is defined in accordance with (4.7) by 
6(ci - di) 



(p,{z) := 



1 



z e 



As in (3.10), the contintuity condition (3.6| at Vi = leads to 
(5.2) ipi{z) — ipi(^ — z^ , for all z e 



and 1 < z < A^. 



Hence, (5.1) and (5.2) imply that the continuity condition (3.6 1 is satisfied if 
ci - di Ci- di 



(5.3) 
and 
(5.4) 



for all 1< « < Af, 



Ici-di Ci - di 1, 1 / • / Ar 

ci\l — Ci\l for all 1 < i < A*. 

aici 



Then, (3.12) implies that the Kirchhoff condition (3.7) is satisfied if 

L N 

(5.5) E^=E^- 

4=1 ' j = L + l ^ 



Similarly as in (4.9), v is solution of (BBMG) with di 7^ if and only if 



di 



(BBMG') 



is a solution of the modified system 

dfUi — aidfdtUi + biUidiUi — on each for t > 0, 
Uj{v,t)^Uk{v,t) fori>0, Vj, k e N{v), 

L,jajdjUj{v, t)=0 for i > 0, V j e Ar(v), 

Thus, without any loss of generality, we can suppose 
(5.6) di = for all 1 < i < n. 

Summing up we have obtained. 
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Theorem 5.1. Let (5.6) hold. If the coefficients > and hi G M* satisfy the 



compatibility conditions 
(5.7) 



— = > for all l<i<n 
ai hi 



(5.8) 



L N 



1=1 



j=L + l 



then there exists a solution u of (BBMG) of the form 

U\^, {t, x) = ip{xi - Cit + Ti), 



where Lp is defined as in (4.7) with 
(5.9) 



Ci > and Ci = 



-Ci. 



Proof. Combining (5.7) and (5.9), we obtain (5.3) and (5.4) with di — for all 



1 < i < N . Then, using the definition of the propagation speeds and (5.7), we have 
Oi _ Oi^/al _ _ ^/oT / hi 



/OiCi 



Cl 



ci 



Thus, up the positive constant (|5.8|) is equivalent to (|5.5|). 



□ 



In the case = 3, recalling the results by Bona and Cascaval in [S], when the initial 
data Wo is the initial profile of a wave 



Wo I,, {xi) = Piixi), for all x^ e e.,, 



for alll < i < iV. 



then the unique solution of (BBMG) is the solitary wave built in Theorem 5.1 



Moreover, our compatibility conditions ^5.71 and (5.8) do not seem to be artificial 



in view of the numerical computations in [8j which exhibit a reflected wave in the 
case iV = 3 and a,- = = di — 1 for all i. 



Remark 6. In view of (5.3) and (5.4), if all the coefficients a^, hi and di are equal, 



i.e., if they all agree with some common value a, b and d, then all the propagation 



speeds are equal. Thus the Kirchhoff condition (5.5) is satisfied if and only if 

N -L = L. 



6. BBM EQUATION ON A TREE 

Now, we consider the case where the graph is a directed tree without boundary 
conditions at boundary vertices. We do not regard our tree as rooted, and in par- 
ticular at each edge there may be more than one incoming edge. 

As in the previous section, we want to construct a solution which is a solitary wave. 
Since paths with more than two edges can occur, we need to add a parameter t; in 
it|^ to account for the edge lengths. Thus, we look for a solution u in the form 

Ui{t,x) ^ ipi{x - Cit + Ti), i = l,...,n. 



where is defined by (4.7) 



From the continuity condition (3.6 ), we deduce that for any vertex Vfc and any edges 



and such that L^i = 1 and i^j — —1, we have 

ipi{l^ - Cit + Ti) = 'fji-Cjt + T,). 
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Figure 5. A finite tree. 



Hence, adjusting (5.3) and(5.4), we need to satisfy the foilowing conditions 

Co 



(6.1) 



^3 



for all i,j e iV(vfc), for all Vfe e Vr. 



(6.2) 

and 

(6.3) 



for all iV(vfc), for all € Vr 



UjCj 



-Ti for all i,j e Af(vfc), for all v^, e Vr 



Up to relabelling we can of course assume that Vq is the root of the tree. Moreover, 
we can choose tq = assuming Vq € eg. Since our graph is a tree, each vertex 
Vi is linked to Vp by exactly one path. Let us denote it by (ej^ , , • • • , e.^^ ) 
for some k € N. Thus, using (6.3), can be uniquely determined by the lengths 
(£ij , ■ ■ ■ , £i^) and by the coefhcients a^^ , bi^ , Ci- and di^ appearing along the path 
from Vq to Vj. As in the previous section, we can without loss of generality assume 
that (5.61 holds. Following the idea of Theorem 5.1 we obtain: 



Theorem 6.1. Let (5.6) hold. If the coefficients > and hi £ M* satisfy the 
compatibility conditions 

(6.4) 



— = > for all i,j G N{\/k), for all Vk G Vr 
a, bi 



and 
(6.5) 



i-kibt = 0, for all € Vr 



then there exists a solution u of (BBMG) of the form 

U\^, {t, X) = Lp{Xi - Cit + Ti) 

where is defined by (4.7) and the propagation speeds are given by 



(6.6) 



Co > and Ci 



—Cj for all Vfe G Vr and all i,j G N{\/k). 



Moreover, the parameters Ti are defined by 



(6.7) To = 



id Ti = 



-Tj + £j for all i,j G N{vk), for all v^, G Vr- 
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Proof. As in the proof of Theorem 5.1 (6.4) and (6.61 imply (6.2) and (6.1) when 
di = 0, and since all the vertices (and edges) are connected to each other, the prop- 



agation speeds are well-defined recursively, starting from cq. Then, (6.7) permits 



to compute all the ti starting from tq since our graph is a tree. (Observe that 
even if the graph is infinite, the path from any given vertex to vq has certainly 
finite length.) Finally, (6.5) is equivalent to the Kirchhoff condition (3.121) when we 



re-write — in the form 



-bj for an arbitrarily chosen iq e N(yk). 



□ 



Remark 7. // the graph is finite and therefore we are in the setting considered 
in [5], then the wave constructed in Theorem 



6.1 is necessarily the unique solution 



of (BBMG), provided that the initial data uq is a solitary wave itself. 



Remark 8. In view of (6.1) and (6.2), if all the coefficients ai, hi and di are 



respectively equal to a, b and d, then all the propagation speeds agree with a common 
value c. Thus the Kirchhoff conditions (6.5 I are satisfied if and only if 



card {i / Lkt = 1} = card {i/ 



-1} 



for all Vfe e 1^. 



In particular, recall this condition is satisfied (in fact, it is equivalent) to the directed 
graph being Eulerian (see e.g. [Ill Thm. 4. 4] J provided that VJ, — 0, i.e., that each 
vertex is of ramification type. 

7. Networks with circuits 

In this section, we consider networks which contain circuits. We first treat the case 
of a graph having one directed circuit, i.e., a path linking a vertex v ^ V to itself 
following the incidence and having more than one edge (see Fig [6]). 




Figure 6. A directed circuit. 



As in Figure [6] we denote by ei, 62, • • ■ , e„ a directed path joining a vertex Vi to 
itself, and let Vi, V2, • • • , v„ the vertices of this path. We look for a solution u of 



(BBMG) in the form Ui{t, x) = ipi{x — Cit + Ti) on each edge e^. According to (3.10 ), 
along the directed circuit, we have 



(7.1) 



^i+liz) = ipi it + T, 



-n+i 



W+l 



for all z £ 



whence 



"jfi+l ( 4+1 + n+l H z 

Ci+2 C.i+2 



I ii H ii+i +n Ti+2 H z 1 for aU z G 

By induction on the length of whole path we can thus prove that 



(7.2) 



tpi{z) — (fi { z + —ti j for all z G 



Thus any travelling wave is necessarily periodic and in particular we obtain: 
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Lemma 7.1. // the graph contains a directed circuit, then there exist no solitary 
wave solutions o/(BBMG). 

Next, we consider the case where the graph contains circuits, but no directed circuits 
(see Fig[7|). 




Figure 7. A circuit which is not a directed circuit. 



It turns out that also on a graph containing undirected circuits a certain com- 
patibihty condition relating the lengths of the different paths between two vertices 
has to be satisfied, in order for a travelling wave to exist. To begin with, we discuss 
the following simple example. 

Example 1. Let us begin by considering the simple case of the graph G in Figure^ 
It is natural to address the following question: After splitting the incoming solitary 
wave in two waves at Vi along the two paths (61,62) and (64,63), can we adjust the 
propagation speeds ( and hence find suitable coefficients of (BBM) ) so that the two 
waves can eventually he glued to form one single outgoing wave at V3 ? 
In order to answer this question affirmatively we need to show that the hypotheses 
in Theorem 6.1 and a compatibility condition stemming from (7.1) can be satisfied 



simultaneously. At vertex •^3, the continuity conditions along the path (61,62) ana 
the path (64,63) and (3.10) imply 



(7.3) 



C2 C4 \ C3 



When we link Vi to itself, from (3.101, we get 

(pi(z) = iy9i(z + /i H — ^^2 -^3 -^4) for all z G M, 

C2 C3 C4 



which is verified since li + ^/2 = ^Is + according (7.3). Linking V2 to itself, 
we have 



Ci C4 C4 



(7.4) 



V2{Z 

Cl 



/ A- ""1/ ""1/ 

'1 H '2 '4 '3 



C2 



C4 



C4 



for all z e 



which is also satisfied because of (7.3). In the same way, linking V3 and V4 to 



themselves respectively, the compatibility condition is verified if ( 7.3 1 is satisfied. 



Using the definition of the propagtion speeds (6.6), we have 

j-i j-i 1 

flL - TT _ TT / _ 
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and we can re-write ( 7.3 ) as 
(7.5) h + , f^k = 



k + 



-k^ 



/ai 



U + 



--k 



We conelude that (7.5) is a necessary condition for the existence of a solitary wave 
on the graph G. 

In the general case, let us consider a graph G with undirected circuits, but without 
any directed circuit. 

Notation 7.2. Let us denote Vgut the set of all vertex v having at least two (directed) 
paths starting at v and going to the same vertex w. For a vertex € Vout? "we have 
at least two directed paths (whose lengths we denote by n and p — n, respectively) 
ending at w € V along which the incidence factors are all equal to 1 . Gluing them, 
we obtain an (undirected) path (e^^, e^^, • • • , e^^^, e^^^^^, • • • , 6;^) with 1 < n < p 
in N such that w — H and = e^^ H e^^. We denote it by [v^, v^J^. 

Along this undirected circuit, we compute a compatibility condition inspired by 



(7.31 to satisfy the transmission conditions at w: 



which can be re- written using conditions (6.6). Thus, we obtain the following. 



Theorem 7.3. Let (5.6) hold. Suppose that the coefficients > and hi G W 
satisfy the compatibility conditions (6.4) and (6.51 for all v^, G and all i,j G 
N(yk)- Then the following assertions hold. 

(1) In order for a travelling wave solution to exist on G, the additional compatibility 
condition 

" 1 ^1 

(7.6) y^^^=li,= -^^k, forv^eVout and all paths [vi.M^Yl 

j=l V h j=n+l V h 

has to be satisfied. 

(2) Conversely, if (7.6) is satisfied, then there exists a solution u o/ (BBMG) of 
the form 

U\^^ (t, x) = ip{Xi - Cit + Ti), 

where ip is defined by (4.7), the propagation speeds are defined by 

(7.7) Co > 



and Ci ~ — -Cj for all G Vr and all i,j G N{\/k) 



and the parameters Ti are defined by 



(7.8) To = 



id Ti ~ 



—Tj + Ij for all Vfe G Vr and all i,j G N(vk). 



Proof. (1) The claim can be proved by induction along the lines of the discussion 
in Example [l]- we omit the details. 

(2) In order to prove the converse implication, we first observe that using conditions 
(6.4), (|6.5[) and (|7.8[) we can construct as in Theorem 6.1 a wave satisfying (BBMG). 



We just have to check that there is no continuity jump in u when considering circuits. 
Let us consider a circuit [v^, v^Jp for some G Vout and let v = f^^ik+i be a vertex 
belonging to this circuit. We want to know whether ipk is well defined when we leave 
V following the paths (ejj._^j, • • • , 6;^) and (e^^, • • • , e^^.). Up to a relabeling we 
can assume that k + 1 < n. As in (7.4 1, from 

'Pij {U, - Ci^ t + Ti. ) = Lpi.^^ {"Ci^^^t + Ti^^^ ) for 2 < j + 1 < n 
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and 

'PiHi~c^J + T,^) = VtH+ii^tH+i - Cz^+it + n^^,) for n + 2<h + l<p , 
we obtain 



\j=k *J "'j 



as well as 



'fc-i 



1 ^2 



Thus, using in particular the compatibility conditions 
and 

that arise from (3.8) we are led to 



fc-i 



yj=fc J 



J = l 



Thanks to (7.6 1 and by definition of the propagation speeds (7.7), this equation is 
satisfied and the wave is well defined along each circuit. The special case k ^ n, i.e. 
V = n Gi^^j , is the condition to glue back the waves into one single wave when 
leaving the circuit. □ 



References 

[1] R. Adami, C. Cacciapuoti, D. Finco, and D. Noja. Stationary states of NLS on star graphs. 
arXiv:1104.3839, 2011. 

[2] V. Adler, B. Giirel, M. Giirses, and I. HabibuUin. Boundary conditions for integrable equa- 
tions. Journal of Physics A: Mathematical and General, 30:3505, 1999. 
[3] H. Amann, Ordinary Differential Equations, de Gruyter, Berlin, 1990. 
[4] J. von Below, Parabolic network equations, Tiibingen, 2"'*edition, 1994. 

[5] J. von Below. Front propagation in diffusion problems on trees. In C. Bandle et al., editor, 
Calculus of variations, applications and computations (Proc. Mont-d-Mousson 1994), volume 
326 of Pitman Res. Notes Math. Ser., pages 254—265, Longman, Harlow, 1995. . 

[6] T.B. Benjamin, J.L. Bona and J.J. Mahony, Model equations for long waves in nonlinear 
dispersive systems, Phil. Trans. R. Soc. Lond. A272 (1972), 47-78. 

[7] H. Berestycki, M. Labadie, and F. Hamel. Generalized traveling waves on complete Riemann- 
ian manifolds, hal-00465620, 2010. 

[8] J.L. Bona and R.C. Cascaval, Nonlinear dispersive waves on trees, Canadian Applied Math- 
ematics Quarterly 16 (2008), 1-17. 

[9] R. Camassa and D.D. Holm. An integrable shallow water equation with peaked solitons. 
Physical Review Letters, 71:1661-1664, 1993. 
[10] S. Cardanobile and D. Mugnolo. Analysis of a Fitzhugh-Nagumo-Rall model of a neuronal 

network. Math. Meth. Appl. Sci., 30:2281-2308, 2007. 
[11] G. Chartrand, L. Lesniak, and P. Zhang. Graphs and digraphs. CRC Press, Boca Raton, FL, 
2010. 

[12] S.N. Chow and J.K. Hale, Method of Bifurcation Theory, Springer- Verlag, New York, 1982. 
[13] A. de Pablo and J.L. Vazquez. Travelling waves and finite propagation in a reaction-diffusion 
equation. J. Differ. Equ., 93:19-61, 1991. 



18 



DELIO MUGNOLO AND JEAN-FRANQOIS RAULT 



[14] M. Ehrnstrom and H. Kalisch. Traveling waves for the Whitham equation. Differ. Int. Equa- 
tions, 22:1193-1210, 2009. 

[15] P.G. Estevez, §. Kuru, J. Negro and L.M. Nieto, travelling wave solutions of the generalized 
Benjamin-Bona-Maiiony equation, Chaos, Solitons and Fractals 40 (2009), 2031-2040. 

[16] J.D. Evans, G.C. Kember, and G. Major. Techniques for obtaining analytical solutions to the 
multicylindcr somatic shunt cable model for passive neurones. Biophys. J., 63:350—365, 1992. 

[17] P. Exner. Momentum operators on graphs. arXiv;1205.5941, 2012. 

[18] J. Guckenheimer, P. Homes, Nonlinear Oscillations, Dynamical Systems and Bifurcations of 

Vector Fields, Springer- Verlag, New York, 1999. 
[19] Y.V. Pokornyi and A.V. Borovskikh. Differential equations on networks (geometric graphs). 

J. Math. Sci., 119:691-718, 2004. 
[20] K. Ruedenberg and C. W. Scherr. Free-Electron Network Model for Conjugated Systems. I. 

Theory. J. Chem. Phys., 21:1565-1581, 1953. 
[21] M. Song and C. Yang, Exact travelling wave solutions of the Zakharov-Kuznetsov-Benjamin- 

Bona-Maliony equation. Applied Mathematics and Computation 216 (2010), 3234-3243. 
[22] G.B. Whitham, Linear and Nonlinear Waves, Wiley, New York, 1974. 

[23] E.T Whittaker and G. Watson, A course of modem analysis, Cambridge University Press, 

Cambridge, 1988. 

Delio Mugnolo, Institut fur Analysis, Univbrsitat Ulm, Helmholtzstrasse 18, 89081 
Ulm, Germany 

E-mail address: delio.mugnolo@uni-ulm.de 

jEAN-FRANgois Rault, LMPA Joseph Liouville ULCO, FR CNRS Math. 2956, Univer- 
siTE Lille Nord de France, 50 rue F. Buisson, B.P. 699, F-62228 Calais Cedex, France 
E-mail address: jfrault91mpa.univ-littoral.fr 



